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We present results of numerical and experimental investigation of the electric breakage of 
a cellular material in pulsed electric fields (PEF). The numerical model simulates the con- 
ductive properties of a cellular material by a two-dimensional array of biological cells. The 
application of an external in the form of the idealised square pulse sequence with a pulse 
duration i^, and a pulse repetition time At is assumed. The simulation model includes the 
known mechanisms of temporal and spatial evolution of the conductive properties of different 
microstructural elements in a tissue. The kinetics of breakage at different values of electric 
field strength E, ti and At was studied in experimental investigation. A 5 mm x 55 mm 
cylindrical slab of apple is taken as a sample. The results of the experimental and numerical 
studies were compared. We propose the hypothesis for the nature of tissue properties evolu- 
tion after PEF treatment and consider this phenomena as a correlated percolation, which is 
governed by two key processes: resealing of cells and moisture transfer processes inside the 
cellular structure. The breakage kinetics was shown to be very sensitive to the repetition 
times At of the PEF treatment. We observed correlated percolation patterns in a case when 
At exceeds the characteristic time of the processes of moisture transfer and random perco- 
lation patterns in other cases. The long-term mode of the pulse repetition times in PEF 
treatment allows us to visualize experimentally the macroscopic percolation channels in the 
sample. We observe considerable differences between the damage kinetics at long and short 
repetition times both for experimental and simulation data. 
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Notation 



A ~ 4dg, cross section area of a single cell or a membrane, m 2 

c = r{ jr l c = <j % J<j{ < 1, resistance moisture transfer coefficient 

C m specific capacity of membrane, F m~ 2 

C* = G m {e w /e m -l)/{2 1 ) 

d m membrane thickness, m 

2d c cell diameter, m 

E electric field strength, kV cm" 1 

E* = 2d c E/u , normalized electric field strength 

AF* = ttuj 2 I (kT-f) , reduced critical free energy of pore formation 

G conductivity of the bond in the network, Ohm -1 

fc Boltzmann constant, 1.381 x 10~ 23 J K _1 

L = 2d c N , total thickness of a sample (slab of apple) , m 

171 = r m/ r m = <T L/ (J rrn membrane resistance resealing coefficient 

n number of pulses 

NxN dimensions of a 2D lattice 

P degree of biological tissue damage 

r resistance of the model resistor in the network, Ohm 

r rn membrane part of r resistance, Ohm 

r c cellular part of r resistance, Ohm 

S(u) survival probability function 

ti pulse duration, /is 

dt impact time duration, or "elementary" time step dt ~ 0.1£^, s 

At pulse repetition time, ms 

T temperature, K 

u transmembrane voltage, V 

u midpoint of a survival probability function S[u), V 

u* = u/u a , normalized transmembrane voltage 

Am width of a survival probability function S(u), V 

Am* = Au/u D , normalized width of a survival probability function S(u) 

U external voltage, V 

W moisture content, % 



Greek letters 



e w = 80, dielectric constant of water 

e m =2, dielectric constant of membrane 

7 surface tension of membrane, N m -1 

A adjustable relaxation parameter 

a conductivity, S m _1 

cr c conductivity of cellular material (barring membrane), S m _1 

<r m conductivity of membrane, S m _1 

r m lifetime of a membrane, s 

Too parameter, lifetime of a membrane at T — oo, s 

T r characteristic time of membrane resealing, s 

Td characteristic time of a moisture transfer processes after PEF treament, s 

lo linear tension of membrane, N 



Superscripts 



2 



d damaged 
i initial 
/ final 

Subscripts 
c cellular 
e effective 
i intact cell 
j juice 
m membrane 
t total 

Abbreviations 

CEF continuous electric field 
PEF pulsed electric field 

I. INTRODUCTION 

Among different nonthermal processing methods used in food technologies, the pulsed electric field (PEF) treatment 
is one of the most promising. A number of new PEF applications were demonstrated for anti-microbial treatment of 
liquid foods, e.g., fruit juices, milk etc., (Barbosa-Canovas, Pothakamury, Palou & Swanson, 1998; Barsotti & Cheftel, 
f998; Wouters & Smelt, 1997), and for the cellular tissue materials (Knorr & Angersbach, 1998, Knorr, Geulen, Grahl 
& Sitzmann, 1994). For years back, the continuous electric field (CEF) treatment was also shown to be good for juice 
yield intensification and for increasing the product quality in juice production (Bazhal & Vorobiev, 2000; McLellan, 
Kime & Lind, 1991; Scheglov, Koval, Fuser, Zargarian, Srimbov, Belik et al., 1988), processing of vegetable and 
plant raw materials (Papchcnko, Bologa & Berzoi, 1988; Grishko, Kozin & Chcbanu, 1991), processing of foodstuffs 
(Miyahara, 1985), winemaking (Kalmykova, 1993), and sugar production (Gulyi, Lebovka, Mank, Kupchik, Bazhal, 
Matvienko et al., 1994; Jemai, 1997). But all these CEF applications were restricted by high and uncontrolled 
increases in food temperature. 

Extension of different PEF applications to nonthermal processing of heterogeneous food materials is limited today 
by the absence of criteria for choosing optimal parameters of PEF treatment and the unclear mechanism of electric 
breakdown processes in the cellular systems. Recently, a significant advance in understanding of the nature and mech- 
anisms of electric field influence on different animal, plant, and microbial cells has occurred (Chang, Chassy, Saunders 
& Sowers, 1992; Weaver & Chizmadzhcv, 1996). The strong electric field causes electroporation of cells, increase 
of their permeability, and, in some cases, disruption of their structural integrity (Zimmermann, 1975). The PEF 
parameters (field strength E, pulse duration ti and number of pulses n) can influence both the degree of membrane 
destruction or structural alteration and the density of pores in membrane (Rols & Teissie, 1998; Gabriel & Teissie, 
1999). The electroporation became very popular because it was found to be an exceptionally practical way of trans- 
ferring drugs, genetic material (e.g. DNA), or other molecules inside the cells (Chang, Chassy, Saunders & Sowers, 
1992; Neumann, Kakorin & Tcensing, 1999). This phenomena is sometimes also called as electropermeabilization. 

For complex material, such as tissue, cellular or food material, PEF application results in increase in the electric 
conductivity and permeability of the whole sample. But the nature of electropermeabilization of complex cellular 
materials is not yet well understood in all details. The cellular materials are highly heterogeneous and the electrical 
properties of such systems depend on the electrical properties of single cells as well as on the geometrical and topological 
properties of materials (Sahimi, 1994). Here the percolation phenomena may play an important role for interpretation 
of the observed experimental results. Particularly, there is no simple relation between the degree of material damage 
and its electrical conductivity. 

A number of ambiguous and yet unexplained phenomena are observed in this field. For example, long-term changes 
in the conductivity of a cellular material after its electric field treatment arc usually observed. It was reported, for 
example, that the conductivity of the vegetable tissue can decrease after termination of electric treatment over at least 
24 hours (Kulshrestha & Sastry, 1998). The explanation of this phenomenon is not trivial since multiple mechanism 
can be responsible for time evolution of the conductivity. 

One of possible explanations is based on the assumption about existence of a partial membrane rcscaling after its 
breakage in a high electric field. Note, that the electric damage of a cell is itself a rather complex process and there 
may exist different time scales for the kinetics of pore evolution. At low strength electric field, E (< 200 V cm -1 ), and 
short pulse duration, ti (~ 10~ 5 — 10~ 6 s) the electrical breakdown is spontaneously reversible, and all the damages 
disappear after the field is switched out (Abidor, Arakelyan, Chernomordik, Chizmadzhev, Pastushcnko & Tarasevich, 
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1979; Weaver & Chizmadzhev, 1996). At moderate PEF treatment (E = 0.5-2 kV cm" 1 , U ~ 1CT 4 - 1CT 5 s) the 
integrity of cells drops rapidly, but due to resealing or recovering process some of the cells loose their permeability 
and the pores may persist in the membrane at larger after PEF application. The resealing process time constant, 
T r may be very large, of order 1 - 10 2 s at 25 °C (Neumann & Boldt, 1990; Chang & Reese, 1990; Chizmadzhev, 
Indenbom, Kuzmin, Galichenko, Weaver & Potts, 1998; Neumann, Toensing, Kakorin, Budde & Frey, 1998; Weawer, 
Pliquett & Vaughan, 1999; DeBruin & Krassowska, 1999). For vegetable food materials the reported resealing time 
constant was of order 1 s (Knorr, Heinz, Angersbach & Lee, 2000). The high PEF treatment (E = 10 — 50 kV cm -1 , 
U ~ 10 -6 s) causes the irreversible damages and this mode of treatment is used for inactivation of microorganisms 
(Barbosa-Canovas et al., 1998). The time constant of a resealing process is a complex function of E and U values 
and depends on the type of cells or membranes. A number of mechanisms were proposed for explanation of resealing 
processes, but in a general case the nature of the long-lived permeabilization is still unclear (Saulis, 1997; Teissie & 
Ramos, 1998; Chizmadzhev et al, 1998; Weaver et al., 1999). 

Another cause of the long-term changes in the conductivity may be related with the different transport phenomena 
in a structured cellular material (Aguilera & Stanley, 1999), e.g. diffusional motion, osmotic flow and redistribution 
of moisture inside the sample, which can be enhanced by PEF application. We can estimate the time constant of the 
diffusion processes inside the cellular material as ~ d 2 c /(QD) w 1 s at 25° C, where d c ~ 10~ 4 m is a radius of cell, 
and D ~ 10 -9 m 2 s _1 (Gekkas, 1992) is an effective diffusion coefficient for the moisture inside a cellular material. 

The aim of this study is to elucidate the mechanism of PEF treatment and long-term changes in the conductivity 
of cellular materials. We consider the damage of a biological tissue in the electric field as a correlated percolation 
phenomena, which is governed by the resealing and moisture transfer processes. The developed simulation model 
includes the known mechanisms of temporal and spatial evolution of the conductive properties of microstructural ele- 
ments (cell membranes, tissue frameworks, etc.). The results of the experimental and numerical studies are compared 
and possible scenarios of the tissue conductive properties evolution after PEF treatment are discussed. 

II. MATERIALS AND EXPERIMENTAL METHODS 
A. Materials 

Freshly harvested apples of the Golden Delicious variety were selected for investigation and stored at 4°C until 
required. A moisture content of apples W was within 80-85%. Typical high resolution scanning electron micrograph 
of the apple sample is presented in Fig. [|. Images were obtained on the instrument XL30 ESEM-FEG (Philips, 
V=15 kV, P=3.5 Torr). The initial specific conductivity of samples (before treatment) was within the interval 
o* = 0.003 — 0.007 S m _1 . The final specific conductivities (after treatment) depend upon the mode of treatment and 
they were within the range of a? = 0.035 — 0.070 S m _1 . The specific conductivity of the apple juice extracted from 
the sample apples was tij = 0.22 ± 0.05 S m . 

B. Experimental methods 

The conductivities were measured by contacting electrode method with an LCR Meter HP 4284A (Hewlett Packard, 
38 mm guarded/guard Electrode-A HP 16451B) for thin apple slice samples at a frequency of 100 Hz and with a 
Conductimetre HI8820N (Hanna Instruments, Portugal) for the apple juice samples at a frequency of 1000 Hz (these 
frequencies were selected as optimal in order to remove the influence of the polarizing effects on electrodes and inside 
the samples). 

Figure ^ is a schematic representation of the experimental pulsed electric field treatment set-up. A high voltage 
pulse generator, 1500V-15A (Service Electroniquc UTC, France) allowed to vary ti within the interval of 10 — 1000 
/is (with precision ±2 /j,s), At within the interval of 1 — 100 ms (with precision ±0.1 ms) and n within the interval of 
1 - 100000. 

Pulse protocols and all the output data (current, voltage, impedance and temperature) were controlled using a data 
logger and a special software HPVEE v. 4. 01 (Hewlett-Packard) adapted by Service Electronique UTC, France. The 
temperature was recorded in the on-line mode by a thermocouple THERMOCOAX type 2 (AB 25 NN, ±0.1°C). 

The thin apple slabs (of thickness 5 ± 0.2 mm and of diameter 55 ± 0.5 mm) were used as samples in the present 
investigation. The freshly cut cells at the outer boundary of a sample cause the initial time dependence of a sample 
conductivity. During the period about 200 s (Fig. ||) the conductivity achieves about 90 % of its stationary value. 
So, in all our experiments we always skipped this transition time before the PEF treatment. Usually the sample 
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conductivity was measured 10 s after finishing of treatment procedure. All experiments were repeated at least five 
times. 



III. DESCRIPTION OF THE SIMULATION MODEL 

A. Probability of a single cell damage 

Weaver & Chismadzhev (1996) gave the comprehensive discussion of different models of the membrane rupture. 
The model that seems to be most reasonable from the physical point of view, is the, so called, transient aqueous pore 
model, in which the average membrane lifetime r m can be estimated with the help of the following expression: 

r m (u) = r 0O exp(AF7(l + U 2 C*)). (1) 

where Too is a parameter (equals to a lifetime of a membrane at infinite temperature, T = oo), AF* = ttoj 2 / (kTj) 
is a reduced critical free energy of pore formation, k = 1.381 x 10~ 23 J K -1 is the Boltzmann constant, a; is a 
linear tension of membrane, N, 7 is a surface tension of membrane, N m , u is a transmembrane voltage, V, 
C* = C m (e w /e m — l)/(2-f), C m is a specific capacity of membrane, F m -2 , e w = 80 is a dielectric constant of water, 
e m = 2 is a dielectric constant of membrane. 

Then the survival probability for a membrane (as a whole cell) during the impact period of dt may be estimated as 

S(u) = exp(-tft/T m (u)). (2) 

Taking Eq. (@) into account, we can rewrite Eq. (|2|) in the following convenient dimensionless form 

S(u*) = exp (- ln2/expa([l - (1 - u* 2 )/(aAu* \n2)}- 1 - 1)) , (3) 

where u* = u/u a , Au* = Au/u a , u Q = (AF* /a - 1)/C*, Au = u Q / ((1 - a/AF*)a In 2), and a = ln(di:/( T oo ln2)) is 
a parameter. 

There is no first principle basis for correct estimations of the different parameters in Eqs. (|l])-([|) (Weaver & 
Chismadzhev, 1996), so the numerical values obtained from fitting of T m (u) to experimental data are used as a rule. 
For example, Lebedeva (1987) presented the following estimations for the lipid membranes: E~ 3.7 x 10~ 7 s , 
w ~ 1.69 x 10~ n N, 7 » 2 x 10~ 3 N m~ 2 , C m ~ 3.5 x 10~ 3 F m~ 2 at 25 °C. ^From these estimations we can obtain 
the following parameters for Eq. (^): u Q ~ 1.52 V, Am* ~ 1.07 and a ~ 1.36 (at dt = 1 /is) and u Q ~ 0.71 V, 
Au* ~ 0.26 V, and a ~ 5.97 (at dt = 100 /is). But for real cellular systems the parameters of Eq. (||) are not clear-cut 
and they would depend on physical properties, type and quality of raw materials as well as on the value of dt. Yet, 
for dcfinitcncss in the following computation we always use the parameters of , u Q and Au* estimated on the basis 
of aforecited data of Lebedeva (1987). 

The example of the survival curve S(u*) is shown in Fig. [|. We see that S(u*) is a kind of probability transition 
function and u* = 1 (u = u a ) corresponds to the midpoint, where S(u) = 1/2. In fact, the value of u may serve 
as an estimate for the critical value of transmembrane voltage, which causes the abrupt decrease of the membrane 
lifetime. Here we define the width of this transition function Am* by drawing a tangent straight line to a curve S(u*) 
in the midpoint u* = 1 as it is shown in the Fig. [| (see dotted line). The dashed line at this Figure corresponds to 
the normalised density distribution function S'/S'max, where 5" = dS/du*. 

B. Simulation procedure 

1. Resistor network model 

We simulate the conductive structure of a cellular material as a two-dimensional array of cells located at the nodes 
of a simple square lattice. The lattice has a size iV 2 with periodic boundary conditions in the x direction in order to 
reduce the finite size scaling effects (Watanabe, 1995). The boundary conditions for y direction are as follows: at y = 
and y = N + 1 we put two electrodes with a constant potential difference U (Fig. |s[) . So the mean drop of potential 
per cell is equal to u = Uj (N + 1) and the reduced voltage on membrane in Eq. (|3[) is defined as u* = U /(u a (N + 1)). 
Then mean strength of the electric field along y-axis is equal to 

E = U/(2d c (N + 1)) =u*u /(2d c ), (4) 
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where 2d c is a cell diameter. 

^From this equation we obtain, for example, E = u*u Q /2d c = 50w* V cm -1 at u a = 1 V and d c = 100 /zm. We 
can introduce also the normalised field strength defined as E* — 2d c E/u Q = u* . As it was mentioned above, the 
exact value of u Q unknown. So we can chose the u Q parameter from the condition of best fitting to the observed 
experimental data. In our simulation we use N — 250 and the total sample thickness is L = 2d c N (ss 5 cm when 
d c « 100 fim). 



2. Micro structural conductive •properties 



We suppose that each node is connected with neigbouring nodes through four conducting resistors, which simulate 
the conductive properties of the cellular media microstructural elements. The resistance of such resistors is determined 
by the two constituent parts 

r = r m + r c , (5) 

which correspond to the membrane (r m ) and cellular (r c ) medium contributions, respectively. Here, membrane 
contribution includes the effective conductive properties of the different membranes in the cellular structure (mainly 
plasmatic and tonoplast membrane). Cellular medium contribution reflects both intra- and extra-cellular conductive 
properties of cellular materials. Intra- cellular contribution includes the effective conductive properties of the cytoplasm 
with its organelles (occupies about 10 % of the cell volume), and the vacuole (about 80 % of the cell volume). Extra- 
cellular contribution includes the apparent conductive properties of the rigid cell wall (occupies about 10 % of the 
total volume and its main structural element is cellulose), of pores and intercellular spaces filled with air etc., (account 
for around 20-25 % of the total volume in apple, see Aguilera & Stanley (1999)). 
We estimate the resistance values of r m and r c in Eq. (0) as 

r m = d m /(a m A), (6) 

and 



r c = d c /(a c A), (7) 

where d m is the thickness of membrane, A is a mean cross-section area of a single cell, a m and a c are the conductivities 
of the membranes and cellular material (barring membranes), respectively. 

At the initial stage of simulation we suppose that all the cells are intact and corresponding resistors in the model 
are equal to 

t*=r i m +r l D . (8) 
In this case the effective conductivity of the whole sample a may be calculated as 

i d m +d c d m + d c cr l m d c /d m 



rA d m /a i m + dc/ai 1 + a l m d c /a l c d n 



(9) 



where we take into account that d m <C de- 
li the potentials in all the nodes are known, u xy , then we can easily determine the transmembrane voltages u 
at all membranes in a system. Consequently, we can determine with the help of Eq. (^) which of membranes will 
destroy after the PEF treatment. The conductivity of these membranes after PEF breakage increases considerably 

Om at -> °°)> and SO C => r fn - °)- 

Figure |§(a) presents the case, when the potential difference in the vertical (y) direction \u x . y — u x ,y+i\ exceeds some 
critical value and, as a result, two cells in the vertical direction are damaged. Note, that at all accounts we always 
observe for this model the simultaneous damage of two cells, because two adjacent cells suffer equal voltage loading. 
In this case, we should make the following interchange of the resistors 

r l ^r d = r{t) = r m (t) + r c (t), (10) 

as it is shown at the Fig. ||(a). The similar case for horizontal (x) direction is presented in Fig. 0(b). Here 
corresponds to resistance of damaged cell, and its time evolution may be found with the help of Eqs. 




G 



3. Resealing processes 



The model accounts for the possibility of temporal electropermeabilization as follows. If any membrane is damaged 
then it begins to reseal immediately, and we suppose, that this resealing results in increasing of the r m as 

r m {t)^ri{l-e-^) (11) 

where rl n = d m j '(a^A), <r( n is a final conductivity of a membrane after the complete resealing, and r r is a time 
constant of resealing process. 

We define the membrane resealing coefficient as 

m = rlJrl = a f m /a l m < 1 (12) 



4- Moisture transfer processes 



The moisture transfer processes at different hierarchical levels, such as diffusional migration, osmotic flow and 
redistribution of moisture inside the sample (Aguilera & Stanley, 1999), enhance as a result of PEF application. The 
new conducting channels arise inside the sample and this causes the temporal decreasing of r c value. We approximate 
this evolution as 

r c (t) = r i c -(r i c -4)(l-e- t ^) (13) 

where r[ — d c / (a[A), a[ is a final conductivity of a cellular material after completion of moisture transfer process in 
the sample, and Td is a time constant of this process. 

For the quantitative description of the moisture transfer processes contribution to the change in cellular material 
conductivity we introduce the resistance moisture transfer coefficient c defined as 

c = 4/r l c = all <o{ < 1 (14) 



5. Finite element analysis 



The simulation of temporal evolution of the system requires a knowledge of the potential distribution in the lattice. 
This distribution can be obtained numerically by solving (Lebovka & Mank, 1992) the discretized version of Laplace's 
equation on a lattice with given boundary potentials. For this purpose we have used the successive relaxation scheme 
(Press et al., 1997). We update the chosen site potential y at the n-th relaxation step according to the following 
equation 



it 



I , f G l U x,y\l + G 2 U x-jy + G 3<, y -1 + G A ul + \ y u _ , , 

'« + A l Gl + G 2 + G 3 + G4 1 ( 5) 



where A is an adjustable relaxation parameter and G\ -j- G4 are the conductivities of the bonds which connect the 
chosen site x,y with all its neighbours (Fig. ^). 

The iteration procedure over all sites in the lattice is continued until the maximum of the relative difference between 
potentials in two successive iterations, u^ y/u^y 1 — 1, converges to a small value 5 (= 10~ 3 ). 



6. Simulated properties and main parameters 

The effective media conductivity a was calculated on the basis of r(x, y) values by applying a highly efficient Frank 
& Lobb (1988) algorithm. The total damage degree P was estimated as the membrane damage degree with the help 
of the following relation 

P =( 1 -^p E r m {x,y)/rC] (16) 
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Note that P = 1 when all cells are damaged (r m (x, y)=rf n = 0) and P = when all cells are intact (r m (x, y) = r % m ). 

We assume the pulse application of external electric field in the form of an idealised square pulse sequence with a 
pulse duration ti, and a pulse repetition time At. In order to increase the accuracy of calculation we introduce the 
"elementary" time step dt which is much smaller then the pulse duration fj. In this work we put dt = O.ltj. 

We use in our calculations the following values of parameters: d m = 5 x 1CP 9 m, d c — 1CP 4 m, a l m = 3 x 10 -7 
S m _1 (Kotnik, Miklavcic & Slivnik, 1998), a[ = 0.1 S m -1 (approximately corresponds to the conductivity of the 
absolutely damaged cellular material), and treat m, c, r r and Td as adjustable variables. 

With this sets of parameters we can adjust the experimentally observed parameters a % and <jf (see Section ( II A )). 
For example, in the case when c = 0.1, we obtain a 1 ~ 3.75 x lO" 3 S m" 1 (Eq. |), and a f ~ 0.100 S m at m = 
(i.e., when resealing is absent) and a* ~ 1.07 x 10~ 2 S m _1 at m = 0.5. 

The example of the simulated kinetics of a breakdown is presented in Fig. || During the period of pulse action wc 
observe destruction of a system and increase of P and a. In the interpulse period the system begins to reseal and we 
can observe the partial decrease of P and a. 



IV. RESULTS AND DISCUSSION 



A. Experimental results 

1. Damage kinetics 

Figure^ presents the examples of the experimental curves of apple slabs relative conductivity ja % versus time t 
dependencies at different values of the electric field strength E and pulse protocols: ti = 1 ms, n = 1 — 15, At = 60 s. 
After each pulse application we observe a rather long-time resealing-like behaviour of <jf / a 1 values during the period 
of order 10 s. So, in each case we measured the equilibrium values of a* /a 1 at time t m ks 10 s after each PEF pulse. 
The results of measurement at two pulse protocols ti = 1 ms, n = 1 — 15, At = 60 s (protocol I), and At = 10 ms 
(protocol II) and different values of the electric field strength E = 200 V cm -1 and E — 500 V cm -1 are presented 
in Fig. [l(]. We see, that there exist significant difference between kinetics of a* /a 1 for these two pulse protocols. As 
we have demonstrated before (Lebovka et al, 2000) the pulse repetition time in the interval of At = 1 — 100 ms does 
not influence the <jf / a 1 versus n dependencies essentially. We have enlarged significantly the pulse repetition time in 
the protocol I (At = 60 s) which resulted in the significant elevation of the conductivity curves to compare with those 
obtained for the protocol II. 



2. Visualization of damages 

Figure |ll| shows the photographs that illustrate the macroscopic structure changes of the apple slabs after PEF 
treatment using protocols 1(a) and 11(b). The dark (brown) spots on the slabs treated using the protocol I seems to 
correspond the formation of the moisture-saturated and more conductive channels in the cellular material, which are 
absent in the case of the protocol II . The visually observed behaviour reveals the different modes of cellular material 
breakage. The only difference between protocols I and II is the pulse repetition time At. The considerable increasing 
of At in the case of protocol II allows us to visualize the existence of certain long-time and large scale moisture transfer 
processes. 



B. Numerical results 



1. Simulation of damage kinetics 

Figure |l2|(a) represents some examples of the simulated <jJ jo 1 kinetic curves after application of n (n = 1 — 15) 
pulses of U = 1 ms duration for two different pulse repetition times At = 60 s (solid lines, protocol I) and At = 10 ms 
(dashed lines, protocol II). We have used in these calculations the following values of parameters: m — 0.1, c = 0.1, 
r r = 10 s, Td — 0.5 s and E* = 0.65. For protocol I we have At > r r ,Td, which means that the resealing and mass 
transfer processes are finished during the time interval At and we observe the typical transition process. For protocol 
II we have At r r , r^, and in this case the system is in transient state during the PEF treatment. For this case we 
observe the increase of the relative conductivity during the time period of order Td resulting from moisture transfer 
processes and subsequent decrease of this value due to resealing processes. 
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Figure 

||(b) depicts the (erf /a l ) tm versus n dependencies for the above-mentioned protocols. The calculated values 
of (erf /er z ) tm were "measured" at time t m after n-th PEF pulse applic ation, a nd this procedure corresponds to the 
real experimental measurement procedure described earlier in Section IV A 1. So, we can consider the Fig . |l2|(b) 



as a simulated analogue of Fig. |T(J. The kinetics of (&f /a 1 ) for the protocol I is represented in Fig. |l2|(b) by a 
dashed-dotted line. It is evident from these data that increase of the pulse repetition time At leads to the el evation 



of (erf /er 1 ) tm versus n dependencies in accordance with the experimental observations as discussed in Section IV A 1 



The examples of simulated breakage patterns for two different pulse repetition times At = 60 s (protocol I) and 
At = 10 ms (protocol II) are shown in Figs. |l3|(a) and (b) respectively. Here, each pattern displays only those cells, 
which were broken after the n-th pulse. We can see that the long repetition time pulse protocol I results in more 
extended and spatially correlated damage patterns (Fig. |l^(a)). Dark clusters of the damaged cells show clearly the 
existence of collective percolation phenomena, which are typical for the electrical breakage of inhomogeneous materials 
(Sahimi, 1994). The short repetition time pulse protocol II results in more rare and uncorrelated damage patterns 
(Fig. 0(b)). 



2. Effects of resealing and moisture transfer processes 

The influence of m and c parameters on the "measured" values of (a* /cr l )t m ~>oo after the application of n = 15 
pulses (Tci = 0.5 s, E* = 0.65) are demonstrated in Figs. 14(a) and (b), respectively. Here, Fig. |l4](a) presents 



(<?f / cr l )t m _ aa versus At dependencies for the case of c = 0.1 and r r — 10 s, and Fig. [u](b) presents the analogue 
dependencies for the case of r r = oo (i.e., when resealing is absent). All these dependencies display characteristic 
dispersion behaviour in the range where At ~ r r , and in all the cases we observe increase of (af /(j% m with At 
increase. 

The results for the simulated kinetics of damage degree P and relative conductivity erf /a 1 for different values of Td 
are presented in Fig. |l5|(a) (r r = 10 s) and Fig. |l5](b) (r r = oo). These calculations were performed for the following 
values of parameters: m = 0.05, c = 0.1 s and E* = 0.65. The data show that resealing influences significantly the 
erf /a" 1 versus t dependencies and this effect is the mostly pronounced at large values of r d . At small r d values (< 5 s), 
the (jf/cr^t) curves show the well pronounced maximum, which is practically absent on P(t) curves (Fig. [T^(a)). The 
interesting finding is that there is no direct proportionality between the damage degree P and relative conductivity 
erf /a 1 . We can observe the obvious decrease of erf /er l even in the case when P(t) practically does not change. This 
behaviour reflects the fact that er* ja % in the percolation phenomena depends not only on damage degree P but also 
on the spatial distribution of the damaged cells over the system (Sahimi, 1994). 



V. DISCUSSION 



The main hypothesis of the present work is that the electric field damage of a biological tissue is a phenomena of 
correlated percolation. This phenomena is governed by the two key processes: resealing of cells and moisture transfer 
inside of cellular structure. In the simulation model we try to use the minimal number of parameters in order to 
imitate only the main feature of this very complex phenomena, which comprises different processes at various micro 
and macro-hierarchical levels of membranes, cells, tissue structure, etc. We have found a considerable difference 
between the damage kinetics at long-term and short-term repetition times At (see Fig. |l0|(a) for experimental data, 
and Fig. |l2|(a) for simulation data). 

In the case of long-term repetition times At > tj a damage process in a system has a correlated character. The 
origin of this behaviour is the following. For heterogenous tissue structure in external electric field the largest gradient 
of field strengths arise near the already damaged cells. So the new cells (after the next PEF pulse) are destroyed 
mainly near the previously damaged cells. In this case the damage processes are spatially correlated and we really 
observe this character of damage distribution at the simulated patterns (see Fig. [13](a)). In the case of short-term 
repetition times At < r<j a damage process in a system has a random character. Cells are also destroyed after each 
PEF pulse, but their damage is latent and is "invisible" for the rest of the system. So, in this case the damage 
processes are spatially random (see Fig. |l3|(a)). 

Unfortunately, we are presently unable to make more strict comparison between the theoretical and experimental 
data as far as we have no precise data for the sets of parameters used in the present model. From the technological 
point of view, it is preferable to use such mode of PEF treatment that allows to achieve more homogeneous damage 
of a cellular material and the maximal degree of damage. These conditions are subject for optimization through 
variation of the PEF treatment mode, but we also need here data on model parameters. 
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We should mention also some restrictions of the present model. This model does not take into account rather 
important details concerning the structure of a cellular material, and particularly, the large scale spatial fluctuations 
of electrophysical properties inside the sample. These fluctuations become experimentally visualized after application 



of PEF with the long-term repetition times At (see Fig. 11(a)) as the large scale percolative channels. For this 
reason we do not discuss here the experimentally observed effects of electric field strength E, because this behaviour 
may be extremely sensitive to the above-mentioned spatial fluctuations. Moreover, our model was developed only for 
two-dimensional systems. It allows us to simulate rather large-scale and realistic systems and avoid the well-know 
finite size scaling effects (Watanabe, 1995). But three-dimensional simulation is, of course, more realistic. All these 
restrictions should be overcome in future models in order to attain more profound understanding of tissue damage 
kinetics. 



VI. CONCLUSION 



The breakage of a biological tissue under the PEF treatment may be described as a correlated percolation phenom- 
ena, which is controlled by two key processes: resealing of the cells and moisture transfer inside of cellular structure. 
The breakage kinetics is very sensitive to the repetition times At of PEF treatment. We observe correlated percolation 
patterns for the case when At exceeds the time of moisture transfer processes and in the other case, when At < r^, 
the random percolation patterns are observed. The long-term mode of pulse repetition times in PEF treatment allows 
us to visualize the macroscopic percolation channels in the sample. 
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FIG. 1. Typical high resolution scanning electron micrograph of the apple sample. The "WET" chamber mode was employed 
that allowed observation of hydrated apple specimens in their natural state. 

FIG. 2. Schematic representation of the experimental set-up used in the study of pulsed electric field treatment of apple 
slices. 



FIG. 3. The example of the initial apple slab conductivity changes caused by the freshly cut boundary regions. The dashed 
line shows the transition time which always is skipped before the PEF treatment. 



FIG. 4. Survival probability for a cell S versus normalised transmembrane voltage u* — u/u calculated using Eq. (|3j) at 
Au* — 0.26, = 3.7 x 10 -7 s, and dt = 100 )j,s. Here, the dotted straight line is a tangent line to a curve S(u*) in the 
midpoint u* = 1 (is shown by the small circle) , and the dashed line corresponds to the normalised density distribution function 
S'/S'max, where S' = dS/du*. 

FIG. 5. The two-dimensional model of the cellular material structure. Each cell is represented by a node with four conducting 
bonds. Here d c is a mean radius of a cell, d m is a membrane thickness, L — 2d c N is the total thickness of a sample, and U is 
an external voltage applied to the sample. 

FIG. 6. This explains the procedure of the resistors interchange when the cells are damaged in vertical (a) or horizontal (b) 
directions. Here resistances n and r,i correspond to the intact and damaged cells and are defined by the Eq. (^) and Eq. (p"o|), 
respectively. 

FIG. 7. The calculation of the chosen site potential u" y according to successive relaxation scheme. For this particular 
case the central intact cell is surrounded by three intact and one damaged cells. The bond conductivities are given by 
Gi = (r(x,y)+r(x,y + l))-i = (2r 1 )" 1 , G 2 = (r(x,y) + r(x- hy))- 1 = (2rT\ G 3 = (r(x, y) + r(x, y - l))" 1 = (2rT\ and 
G4 = (r(x,y) + r(x + l,y)) 1 = (r ! + r d ) 1 , respectively. 

FIG. 8. The example of the simulated breakdown kinetics: degree of breakdown P and effective conductivity a of the 
system versus number of pulses n. Here t i = 1 ms is a pulse duration, At = 5tj = 5 x 10~ 3 is a pulse repetition time, r r = 10 -2 
s is a resealing time, and Td is a mass transfer process time. The calculations are performed at E* = 0.75, m — 1 (a case of 
complete resealing) and c = 0.1. Here, all the parameters are chosen only with the purpose of clear illustration of the work of 
algorithm. 

FIG. 9. Examples of relative conductivity or* /(J 1 of apple slabs versus time t dependencies at different values of the electric 
field strength E = 500 V cm' 1 and E — 200 V cm -1 and pulse protocols: ti = 1 ms, n = 1 — 15, At = 60 s. 



FIG. 10. Relative conductivity /a 1 of apple slabs versus number of pulses n at different values of the electric field strength 
E — 200 V cm -1 and E = 500 V cm" 1 and two pulse protocols: ti = 1 ms, n = 1 — 15, At = 60 s (protocol I), and At = 10 
ms (protocol II). In all the cases the value of 0* ja % was measured at time t m = 10s after the end of PEF treatment. All the 
experiments were repeated five times. The error bars represent standard data deviations. 



FIG. 11. Photographs which illustrate the structure changes of the apple slabs after PEF treatment at E = 500 V cm 
ti — 1 ms, n = 10 and different pulse repetition time At = 60 s (a, protocol I), and At = 10 ms (b, protocol II). 
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FIG. 12. Calculated curves of relative conductivity a 1 /a 1 versus time t (a) and number of pulses n (b) for two different 
pulse repetition times At = 60 s (protocol I) and At = 10 ms (protocol II). The following values parameters were used in these 
calculations: U = 1 ms, m = 0.1, c = 0.1, r r = 10 s, Td = 0.5 s, n = 1 — 15 and E* = 0.65. 

FIG. 13. Simulated breakage patterns for two different pulse repetition times At — 60 s (a, protocol I) and At — 10 ms (b, 
protocol II). We have used the following values of parameters in these calculations: U = 1, m = 0.1, c = 0.1, ry = 10 s, r d = 0.5 
s and E* — 0.65. In each pattern we display only those cells, which were broken after the n-th pulse. 

FIG. 14. Calculated (a f /cr i )t m — ► oo versus At dependencies at different values of m (a, C = 0.1 and r r = 10 s) and c (b, 
r r = co, i.e., when the resealing is absent). All calculation were performed at Td = 0.5 s. The relevant repetition time when 
At ~ r r is shown by dashed line. 

FIG. 15. The simulated kinetics of damage degree P and relative conductivity cr^/cr 1 for different values of Td, T r = 10 s 
(a) and r r = oo (b). Arrows show the direction of Td increase. All the calculations were performed for the following values of 
parameters m = 0.05, c = 0.1 s and E* = 0.65. 
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